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Abstract. This contribution briefly illustrates preliminary calculations of the electromagnetic form 



■ factors of ^He and ^H, obtained within the Light-front Relativistic Hamiltonian Dynamics, adopting 
i) a Poincare covariant current operator, without dynamical two-body currents, and ii) realistic 
nuclear bound states with S, P and D waves. The kinematical region of few (GeV/c)^, relevant for 
forthcoming TJLAB experiments, has been investigated, obtaining possible signatures of relativistic 

■ ^ , effects for > 2.5 {GeV /cf. 
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■ INTRODUCTION 

a^ 

^ ! Relativistic Hamiltonian Dynamics (RHD), introduced by Dirac in a seminal paper 

tJ" I represent a viable tool for fulfilling the Poincare covariance, i.e. a fundamental 

■ symmetry. In particular, in the study of the electromagnetic (em) form factors (ff) of 
O I few-nucleon systems some advancements were made within the so-called Light-front 

I (LF) RHD (one out the three ones proposed by Dirac), both calculating the Deuteron em 

observables [2] and obtaining a first description of the Trinucleon ff [3] (retaining only 
S and S' waves). It should be pointed out that relativistic calculations have a twofold 
rN I interest, both from a general point of view and from a phenomenological one, given the 

■ forthcoming experiments at TJLAB, in the region of few {GeV /c)^ (see, e.g., Ref. [0] 
forthe^Heand^Hff). 

Our aim is to construct, within the LF RHD, a relativistic approach for light nuclei, 
taking into account only a fixed number of degrees of freedom, that i) automatically em- 
beds the whole successful phenomenology already developed within the non relativistic 
framework, and ii) includes non perturbatively the relativistic features requested by the 
Poincare covariance. The Bakamjian-Thomas (BT) procedure (see, e.g., Ref. [15Q) allows 
us to implement the above program, since within such an approach it is possible to ex- 
plicitly construct operators that fulfill the commutation rules of the Poincare generators, 
in terms of i) operators that do not contain the interaction and ii) an operator that acts as 
the mass of the interacting system. This interacting mass operator, that depends upon in- 
trinsic variables, has to satisfy the same constraints given by the Galilean group, namely 
the same properties implemented in the non relativistic quantum mechanics. In sum- 
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mary, one can exploit realistic wave functions for light nuclei (for A=3, see, e.g., Ref. 
[0]), that depend upon suitable Jacobi coordinates, in order to evaluate matrix elements 
of a Poincare covariant current operator [7] for an interacting system. The appeal of the 
BT procedure, within the LF framework, is given by the fact that relativistic effects im- 
posed by Poincare covariance can be straightforwardly taken into account through the 
relativistic kinematics and the presence of the so-called Wigner-Melosh rotations (see, 
e.g., Ref. lUl]), that ultimately lead to using the standard Clebsh-Gordan machinery for 
obtaining many-nucleon wave functions with the correct angular coupling. 

FORMALISM 

Among the main features of LF RHD (for an extended review see, e.g., Ref. one 
has to mention the largest number of kinematical Poincare generators (seven, as dictated 
by the symmetry of the initial hypersurface jc+ = 0) and the simplest procedure for 
separating out the center of mass motion from the intrinsic one, in strict analogy with the 
non relativistic procedure (given the absence of the square root in the operator generating 
LF-time translation of the system). On one side, since the LF boosts form a subgroup of 
the kinematical set, it is not necessary any Wigner function for taking care of the boost 
transformation of the nuclear wave function, when we consider, e.g., a nuclear final 
state recoiling in the laboratory frame. On the other side, since the rotations around the 
X- and j-axes are dynamical ones, we need to overcome a difficulty. A possible strategy 
is suggested by the BT construction, that amounts to put in relation the LF spin with the 
canonical one through unitary transformations, namely Wigner-Melosh rotations. For 
the sake of concreteness, the Wigner-Melosh rotations in the 2x2 representation, are 
given by 

= xl'^^^^^^^^^x. (1) 

y'(m + ^,+)' + |k,-x|2 

where L^^{k^) and LLf{k^) are the SL(2, C) representations of the canonical boost and 
the LF one, respectively (see, e.g. Ref. [0]), and is the intrinsic momentum of the i-th 
constituent. 

Therefore, the canonical spin of a single constituent (i.e., the total angular momentum 
in its intrinsic frame) is given in terms of the LF one, by 

sc{i) = [i^w{kj)] slf{i) (2) 

and the component of a state 1 1//) in the LF-spin basis is related to its canonical counter- 
part, by 

LF{o\w)=Y,D^[Ml{k])U, ,(ct'|V/) (3) 
a' 

This illustrates the source of a very cumbersome algebra necessary to perform Trinu- 
cleon calculations. 
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Following Ref. fy*], for an interacting system, an em current operator, 7^(0), that 
fulfills the extended Poincare covariance (i.e. considering parity and time reversal, as 
well) and Hermiticity, can be constructed by a suitable auxiliary operator, j^, that fulfills 
rotational covariance around the z-axis in a Breit frame (P/ + P/ = 0), where the ± 
component of the momentum transfer is vanishing (q^ = 0). Note that such a frame 
is dijferent from the Drell-Yan one, where = 0, where the kinematical symmetry 
around ^ is not exploited. In general, the matrix elements {Pf\J^\Pi), still acting on 
internal variables, are directly given by the matrix elements of the auxiliary operator 
j^, evaluated in the chosen Breit frame. A possible Ansatz for a many-body auxiliary 
operator is built from i) the free current (a one-body operator) and ii) the ± component 
of the angular momentum operator S (a many-body operator in LF) as follows 



(4) 



where rx{Q) is a rotation by an angle 6 around the x-axis, J^^iiqe^ = Tlf Jf^eei^) 
with n the projector onto the states of the (initial or final) system and S = the LF-spin 



operator of the system as whole: it acts on the "internal" space. The operator, Jf^-g^iO) is 



the proper sum over A=2,3... free Nucleon current given by 

4 - -F2'{A') ^^^^^ + f [F^{^') +i^2^(A^)] 



(5) 



where = {p'^ — P^Y and F^2) {^) the Dirac (Pauli) Nucleon ff. It should be pointed 
out that the Nucleon ff depend upon p'^ — p^ and not upon q^ ^ p'^ — p^, since only 
three components of the four-momentum, i.e. /?+ and p^, are conserved quantities. On 
the other hand, within RHD framework all the particles are on their mass-shell. 

In the chosen Breit frame, charge normalization and current conservation (for Mf = 
Mi) can be fulfilled by imposing ^ - (qe,) = ^ + (qe,) tA\^. 

Summarizing, if 7^(0), that is an operator acting in the whole space, is Poincare 
covariant, then the intrinsic operator, j^, is invariant for rotations around the z-axis, and 
viceversa. Moreover, in Eq. dH), Sx generates many-body contributions to j^, as well 
as ^^^{qez), if a many -body term is added to Jf^eei^)' discussed in the following 
Section. 

For evaluating matrix elements of i^{qe^, the eigenstates of the interacting system 
are needed. To this end, one can use the "non relativistic solutions", but with Wigner- 
Melosh rotations in the angular part, if the interaction V = M^t — Mq (where Mint 
is the mass operator of the interacting system and Mq the corresponding free mass) 
can be embedded in a BT framework. The BT construction for obtaining interacting 
Poincare generators suggests a necessary (not sufficient) condition on the interaction 
(see Ref. dsj]): V must depend upon intrinsic variables combined in scalar products, 
i.e. [Mlf:V] = [So,V] = [P^,V] = [P^,V] = 0, where Mlf are the LF boosts, Sq = the 
angular momentum operator for the non interacting case (since Sq = 5? ^ and So,^ = Sim^^^ 
the eigenvalues of Sq and 5o.7 can label the eigenstates of the interacting system). The 
non relativistic interaction fulfills the above requirements. 
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A'/ - k[ Kj - ki K, - k\ K, - ki Kf - k\ K, - k, Kf - k'l K, - fci 

(a) (1)) (O (i) 

FIGURE 1. A class of two-body (at least three particles in flight) current contributions generated 
within a LF analysis of a Yukawa model, see text. Diagram (b) represents a pair production contribution. 
Diagrams (c) and (i) contain the instantaneous, in LF time, propagation of an internal fermionic particle, 
marked by a horizontal dash. The LF time flows from the right to the left. (After ^) 



TRINUCLEON EM OBSERVABLES 

The choice of a Breit frame where = 0, namely 7^ 0, is a far reaching one, since, as 
above mentioned, one can find a simple constraint to be fulfilled by a one-body intrinsic 
current operator for recovering the Poincare covariance, as well as by each many-body 
term that one could add to J^^g^. Furthermore, this choice necessarily produces a two- 
body current related to a pair production [jSO (see diagram (b) in Fig. [T]), compelling us to 
consider the inclusion of a larger set of two-body currents in the future calculations. In 
particular, a recent analysis of a 4D Yukawa model [8], in ladder approximation, has led 
to a 3D current on the LF, fulfilling the Ward-Takashi Identity (for a general discussion 
see itsl]). In Fig. ([T]), a set of contributions to the first-order (in the interaction) 3D 
LF current is shown. In the preliminary calculations of the Trinucleon em observables 
presented in this contribution, the two-body terms, like the ones depicted in Fig. ([I]), are 
not included, while the application to the Deuteron case is in progress. One can anticipate 
that i) the pair term affects all the three Deuteron ff, while instantaneous terms (present 
only for fermionic constituents) contribute to the magnetic one, ii) the pair term vanishes 
for 0, as it must do, while instantaneous ones survive, iii) the remaining, on-mass 
shell terms (like (a) in Fig. ([T])) affect all the ff in the whole range of q^ , iv) the pair 
term should be maximal at ^+ ~ m^/ (cf the discussion in [9]). The macroscopic current 
of the Trinucleon is given in terms of charge and magnetic ff by 

where = ±1/2 labels ^He and respectively. Microscopic evaluations of the ff are 
based on proper traces of the current in Eq. (H)), viz 

= - = \ Tr[tm = 1 rr[^+(r,)] 

F^asiQ^) = f[^~+F^^ = -i^ Tr[dy UTz)] = Tr[ay J^,(rj] (7) 

where the matrix elements of the microscopic current are given by 

^^,^(rj^(vi't;;,,/"|^n^?^^) (8) 



4 



with 1^1 the Trinucleon wave function. In the actual calculation we have used the 
bound states obtained through a variational technique @] by taking into account two- 
nucleon forces, like AVIS [10], and three-body ones, like UIX fll]. Moreover, ^He and 
are distinct, since the Coulomb forces are included. The bound states contain 5, S' , 
P and D waves, and the Wigner-Melosh rotations suitable for expressing the LF spins 
in terms of the standard ones. It should be pointed out that the numerical calculations 
involve 6D Montecarlo integrations. 

In Tables [H and [21 preliminary results for the static em properties of ^He and -^H are 
shown. In particular in the first Table, the bound states obtained by retaining only two- 
nucleon forces are presented, while in the second Table the bound states correspond to 
AVIS + UIX. The benefits of considering the Poincare covariance are clear and of the 
same order found in the Deuteron case [2]. We expect a possible improvement of the 
magnetic moments by including the instantaneous contributions (cf Fig. [T]), since they 
particularly affect the magnetic ff only (as already seen for the Nucleon case [112D). 



TABLE 1. Preliminary calculations of magnetic moments and charge 
radii of ^He and ^H. The two-body force, AVIS, and the Coulomb in- 
teraction are included. Trinucleon wave functions from |j6(]. Probability 
of the waves considered: ^s+s'iAvl8) - 91.4%, ^/.(Avl8) ~ 0.07%, 
^d(Av18) ~ 8.5% 



Theory 


M('He) M('H) 


r,,(%e)fm 


r,/,(3H)fm 


NR(Sh-S') 
LF(Sh-S') 
NR(Sh-S'h-Ph-D) 
LF(Sh-S'h-Ph-D) 


-1.700(1) 
-1.758(1) 
-1.762(1) 
-1.834(2) 


2.515(3) 
2.600(3) 
2.579(2) 
2.674(2) 


1.926(3) 
1.949(3) 
1.916(4) 
1.941(4) 


1.726(3) 
1.771(3) 
1.718(4) 
1.759(4) 


Exp. II -2.1276 


2.9789 


1.959(30) 


1.755(86) 



TABLE 2. Preliminary calculations of magnetic moments and charge 
radii of ^He and ^H. Two- and three-body forces, AV1B: + UIX, are 
included, as well as the Coulomb interaction. Trinucleon wave func- 
tions from [6]. Probability of the waves considered: ^5_|_y(Avl8 + 
UIX) ~ 90.5% ^p{AvlS + UIX) ~ 0.01% 3^d{Av1S + UIX) ~ 9.3%. 



Theory 


Al(^He) 




rc/,(^He)fm 




NR(Sh-S') 
LF (Sh-S') 
NR(Sh-S'h-Ph-D) 
LF (Sh-S'h-Ph-D) 


-1.697(1) 
-1.759(2) 
-1.760(1) 
-1.837(2) 


2.494(2) 
2.588(2) 
2.569(2) 
2.669(2) 


1.848(3) 
1.870(3) 
1.841(4) 
1.867(4) 


1.695(3) 
1.712(3) 
1.666(4) 
1.690(4) 


Exp. 


-2.1276 


2.9789 


1.959(30) 


1.755(86) 



In the evaluation of the Trinucleon ff two different sets of Nucleon em ff have been 
considered: i) the Gari-Kriimpelmann Nucleon ff [13] and ii) the ones obtained within 
a novel LF approach see Figs. [2] and [3l in order to test the dependence of the 
Trinucleon ff upon the new features of the Nucleon ff, like a possible zero in the proton 
charge ff|ll2Lll4|]. 

The charge and magnetic ff of -^H and -^He, evaluated in the Breit frame where = 
and with S, S', P and D waves, are shown in Figs. |4] and [5] for bound states correspond- 
ing to two-nucleon forces, while in Figs. [6] and |7] the calculations with three-body forces 
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FIGURE 3. The same as in Fig.lH but for the magnetic Nucleon ff. 



are presented. Calculations with only 5 + 5' waves are shown in f^. Some expected fea- 
tures, like the presence of large relativistic effects on the tails for 2 > 7 (l//m), are well 
confirmed, as well as the necessity of two-body dynamical corrections to the current op- 
erator (cf the positions of the minima). An interesting signature of the three-body forces 
can be found in the tails, since they give more binding and smaller charge radii. It could 
be relevant, if such an a effect will still be present after including two-body dynamical 
currents. As a final remark, we should note that differences between relativistic calcula- 
tions obtained by using the Gari-Kriimpelmann ff and the ones calculated adopting the 
LF ff are not sizable. 



CONCLUSIONS & PERSPECTIVES 

In order to embed the Poincare covariance in the description of light nuclei we adopt 
a Light-Front RHD and the Bakamjian-Thomas procedure. Extending our approach, 
already applied to the Deuteron case [|2|], the em observables of -^He and have been 
calculated for the first time with all the waves, 5, S', P and D, in the bound states 
and taking into account the three-body forces as well. The relativistic effects on em 
observables at = 0, though of the order of few % but in the correct direction, are 
encouraging. Moreover, the sizable effects at high indicate the essential role played 
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FIGURE 4. Preliminary calculations of and ^He charge ff vs Q^. The adopted Trinucleon wave 
functions ["6^ contain both the AVIS two-nucleon interaction and the Coulomb interaction. 5 + 5'' + P + Z) 
waves are taken into account. Thick lines: LF calculations. Solid line: full calculation and LF Nucleon 
ff 1IT2II : dashed line: full calculation and Gari-Kriimpeknann Nucleon ff lfl3il . Thin lines: non relativistic 
calculations. Data from ifTsll . 




Q (1/fm) Q (1/fm) 

FIGURE 5. The same as in Fig. |4] but for the magnetic ff. 



by the Poincare covariance for analyzing the em ff in the region of few GeV's. Notably, 
three-body forces could be important in the same kinematical region. 

A full calculation, with a systematic inclusion of two-body currents, like the ones 
shown in Fig. [H will be presented elsewhere. 
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